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I . 



INTRODUCTION 



The need for readily computed parameter estimates is 
great. Maximum likelihood estimators are known to be asymp- 
totically efficient, but in many settings they are hard to 
find. The most popular alternative is the method of moments 
which usually yields readily computed estimates, but some- 
times these estimates are not very efficient. See ref. [2]. 
This study looks at the efficiency of a method which uses 
the probability generating function (empirical La Place trans 
form evaluated at a desirable value of its argument, u q ) . 

The method presented herein requires computing power 
greater than that of the method of moments, but less than 
that of maximum likelihood, which calls for a psi function 
capability . It can be used with a hand held calculator that 
has the lower order transcendental functions, i.e., logarithm 
roots and power. 

The basic idea is to select a system of estimating equa- 
tions which equates various statistics to their expected 
values. In chapter II, estimating equations were set and 
efficiency computations were done using theoretical work 
presented in [3] . Chapter III contains the procedure to be 
followed, applications of the method to some data sets, and 
comparisons with the other methods. 

In chapter IV, efficiency and optimum u q and t (=-log u q ) 
tables for various p-r combinations are given and efficiency 
contours in three different planes are graphed. 
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II. METHODOLOGY 



A negative binomial random variable X has the probability 
function 



f (X=x; r , p) = q X p r (2.1) 

for x= 0,1,2, ... 0 < r, 0 < p < 1, p + q = 1. The partial 
derivatives of its logarithm are, using h = log f 

9h _ r _ x 
3p P q 

= V (r + x) - V (r) + log p 

Using the basic recursive formula for the psi function [1] 

x 

Y(r + x) - 'F(r) = l - _ 1 1 + j 

j=l 

one may then express the system of maximum likelihood 
equations as 

x - = o 

p 

(2.3) 

l f^t+t + log p = 0 

n J 



Ave 
i=l / 
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where the values x, , . . . , x are the data that result when 

1 n 

the population (2.1) is sampled n times. The system (2.3) 
is nonlinear in r,p and difficult to solve. Iterative 
methods must be used and the second member of (2.3) (or 
some version thereof) must be recomputed in each cycle. 
This is the main computational difficulty in using maximum 
likelihood in this setting. 

The information matrix A, which is defined as: 



A 





a 2 h 



3r 3p 



3 2 h 

3p9r 



3 2 h 

3 r 2 



is needed for efficiency calculations and can be developed 
using the methods presented in [2] . The result is 



A 



where 



r 



2 

qp 



i 

p 



i 

p 




A 22 - V ' (r) - E(T (r + x) ) 
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The determinant of A is found (using [2]) as 



|A| 



OO 

1 r q 11 r ! n ! 

2 l (n+1) ' (r+n) I 

q n=l 



(2.4) 



The two estimating equations g^ and g^ f which are 
exploited herein, are given as 



g-L (r , p;x) = x - ^ = 0 



g 2 (r,p;x) = t X - ( 1 J^ t ) r = 0 



(2.5) 



where 



i=l 

is the empirical generating function. (Note: the substitution 

t = exp(-u) converts this to La Place transform). Reduction 
to a single equation is effected by solving gj^ = 0 for p 
in terms of x and r. 



p = — (2.6) 

x + r 

and substituting this into g 2> A function f(r) is obtained 
in terms of known variables, namely t , x, and t. 
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f (r ) = ( ) r - t X = 0 (2.7) 

r + x(l - t) 

We can solve equation (2.7) for r by Newton's method, the 
iterative relationship being 

f (r.) 

r i+l = r i - nFTy 

i 



where 



f'(r) = ( 3 ) r • [ ( — ) + log ( ^ )] 

r + x(l-t) r + x(l-t) r + x(l-t) 



assuming an initial value of r and a suitable t can be 
found (see table 4.2). After solving for r, p can be found 
using ( 2 . 6 ) . 

Efficiency for this r-p pair is given by (see [2,3]) 



where 



Eff 



I m " 1 1 
|A| 





( 2 . 8 ) 



(2.9) 



| A | is given in (2.4), the matices A and C in (2.9) are 
defined as 
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A = E 



. 9g: 



3p 



3g, 

3p" 



39j 

3F 



3g, 

3r~ 



C = nE 



<3i<3 



1*2 



g l g 2 



and g x ,g 2 are from (2-5). After some calculation, A and 
C turn out to be 



q 

P 



A = 



( 2 . 10 ) 



-( 



1-qt 



r-1 r(l-t) 



(1-qt) 



p .r . ,1-qt. 

I^t) l0g( p ] 



rq 

2 

P 



_ rqp r 1 (l-t) 



(1-qt) 



r+1 



C = 



( 2 . 11 ) 



rqp r ~ 1 (l-t) 



(1-qt) 



r+1 



1-qt 



r) r " ( 



E_) 2r 



1-qt 



The efficiency of the estimation scheme is found, substituting 
|A| , |c| and j A j into (2.9) and (2.8), as 



12 



Ef £ = 



_ rp 



r+2 



[<A2£> legist) - (2ii^-)] 2 

ir it 



, 2 (r+1) 



q [ (l-qt) “ ^ ’ ~ u; _r ,, r ,, , v 2 n 



n 



(l-qt 2 ) r 



■ - p r (1-qt) ^ - rqp r (1-t) j [ £ 



n=l (n+1) 



(r+nTT 

( 2 . 12 ) 



■] 



Detailed computations are given in Appendix A. The initial 
value of r is found using the method of moments (p,r) , 
and the value of t used maximizes the efficiency (2.12) 
when evaluated at p, r. 
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III. USAGE AND APPLICATION 



In this chapter the procedure to be followed, in order 
to use the method developed, is given in greater detail and 
is applied to three sets of data taken from [4] and sixty 
sets of data generated by computer simulation using variqus 
combinations of the parameters. 

A. PROCEDURE 

The steps of procedure are given below. 

1. Using the method of moments, find p and r as 
starting values. Starting values are given by (3.1) 

(see [2]) 



P 



x/s 2 



(3.1) 




— . 2 . 
where x is the sample mean and s is the sample variance. 

2. Using p, r found in step 1, find the value of t, 
which maximizes efficiency function (2.12), and call it 
t , and also find the efficiency, which is maximum with this 
pair p, r and t , and call it EFF1 . Tables (4.1) and (4.2) 
have been prepared for this. Single variable search and 
golden section search were used by the author to find t Q 
and EFF1. 
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3. Using r as a starting value and t , found in step 2, 
find r* as a solution of (2.7), and p* using (2.6). Newton's 
method was used by the author to find r*. 

4. Step 2 can be repeated with p* and r* replacing p 
and r, and using a new t Q (update t from Table (4.2) if 
the current estimated efficiency (Table 4.1) is not suffi- 
ciently high) . This step involves the recomputations of 

X 

t and is seldom needed. The new efficiency using p* , r* 
and updated t is called EFF2 . 

5. Stop if all the following conditions are met. 

(1 - EFF1 ) < e x , or |EFF2 - EFFl | < e 2 



and 



|P* " P| 1 e 3 / 



and 



l r * " r l i e 4 • 

Then p* and r* are the estimators of true parameters p and 
r, and the estimated efficiency is EFF2. 

6. If one or more of the conditions in Step 5 are not 
met, let 



P 



p* 
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EFF1 



EFF2 



and go to step 3 and repeat the steps following until the 
conditions in step 5 are met. e^, e^ and e^ are the 

stopping criteria which can be chosen by the user . 

B. EXAMPLES 

Examples one through three are the applications of the 
method to the Cricket score data of Reep, Polard and 
Benjamin [4] . Example four is based on data generated by 
computer and contains sixty cases . 

Example 1 : Applying method developed to the Cowdrey 

data [4] the following information is obtained. 

Sample mean = 1.692 

Sample variance = 4.343 

p* = 0.325 

r* = 0.816 

Optimum t = 0.459 

Efficiency = 0.996 

Number of Iterations = 3 

. ... , >\ /% 
With the same data maximum likelihood estimators p, r, 

and method of moment estimates p, r are given in [2] as 



0.329 



0.831 



P 



A 

r = 



p = 0.390 



r = 1.080 



Example 2 ; Example 1 was repeated for Barrington data 
[4] and following information is obtained. 

Sample mean = 2.095 

Sample variance = 4.939 

p* = 0.346 

r* = 1.111 



Optimum t = 0.538 

Efficiency = 0.995 

Number of Iterations = 3 

p = 0.345 r 

p = 0.424 r 



1.014 

1.543 



Example 3 ; Like the previous two examples, the following 
information is obtained by applying the method to the 
Graveney data [4] : 

Sample mean = 1.570 

Sample variance = 4.474 

p* = 0.315 

r* = 0.722 



17 



Optimum t 



o 



0.430 



Efficiency = 0.996 



Number of iterations 



2 



P 



0.317 



r 



0.729 



P 



0.351 



r 



0.849 



When the first three examples are examined carefully 
it is noticed that the method developed in this study is 
almost as efficient as maximum likelihood method. Estimators 
are pretty close to maximum likelihood estimators and 
much better than those found by the method of moments. 

Example 4 ; This example is based on the data generated 
by computer. Sixty cases, each of which has a different 
sample size and parameters p and r. Four sample sizes, 
which are 15, 30, 50 and 100; three r values, 0.5, 2.5 and 
5; and five p values, 0.05, 0.1, 0.3, 0.5, 0.8 are used. 

In each case the data is generated for given sample size 
(n) , p and r. Results, which contain case number (case no.), 
sample size (n) , p, r, method of moments estimates (p,r), 
optimum u Q at p-r (Initial u q ) , efficiency at the end of 
the first iteration (Initial eff.), empirical La Place 
transform estimates (p*,r*), optimum u Q at p*-r* (Final u Q ) , 
efficiency at p*-r*-final u Q (final eff) and number of 
iterations (No. of itr.), are tabulated in the following 
pages. Data generated for the sixty cases are given in 
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Appendix B. The computer program, written in FORTRAN IV, 
is also given in Appendix C. 

When all the cases are examined, it's seen that esti- 
mates are not given for cases 2, 3 and 18. Case 3 provides 

nothing to work with as all samples are zero, yielding 

— 2 2 — 

X = 0, s =0. The other two cases have s < x which 

signals trouble because the method of moments estimates 

r < 0, p > 1; values which cannot be used to initialize 

our iteration scheme (3.1) or the maximum likelihood scheme 

(2.3) . 

For each of these cases the function f (r) of (2.7) 
is a non-negative decreasing function of r which tends to 
+«> as r ■+ 0 and is asymptotically zero as r -*■ 00 . Thus no 
root exists. The maximum likelihood approach yields a 
comparable situation. 



19 



4-1 

o 



.u 

g H 



ro 



ro 



ro 



ro 



ro 



ro 



ro 



ro 



ro 





o 




00 




00 


00 


in 




cr» 


ro 


rH 


00 


ro 


4-1 


o 


, cn 




<J\ 


<J\ 


<J\ 


a> 


00 


00 


00 




00 


r- 


w 


• 




<J\ 




<J\ 


<J\ 


<r> 


cr> 


CP» 


ON 






<j\ 


w 


H 


• 

























CN 


r- 


<3< 


rH 




C7* 


rH 


CN 


ro 


rH 


r- 


CN 


a\ 


^D 






O 


rH 


rH 




CTi 


r- 


CN 


in 




CN 


rH 


i— i 




00 


CN 




00 




o 




CN 




a\ 




CVJ 


CN 


00 


rH 


rH 


in 


o 


rH 


in 


o 


o 





■K 

!h 



•K 

a 



vd 


in 




rH 


o 


CN 


<D 


rH 


cj\ 


in 


ro 


rH 


ro 




CN 


CN 




00 


O 


CN 


CN 


r- 


00 


rH 


in 


CN 




i ^ 


in 


<y\ 


in 


in 


rH 


r- 


o 


<D 


00 


r* 


in 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


VO 1 


1 


ro 


o 


in 




i H 


ro 


ro 


o 




ro 


o 


rH 


o 


a\ 


o 




<J\ 


<j\ 




in 


rH 


x CN 


r^ 


o 


00 


on 


a\ 




CN 


CN 


rH 


r- 


ro 


rH 


in 




in 








in 


ro 




ro 


o 


rH 


rH 


o 


o 


o 



3 



a 



a\ 


00 


00 




00 


00 


rH 


00 


<J\ 


ro 


rH 


00 








G\ 


<j\ 




G\ 


a\ 


00 


00 


00 




00 


VD 


°! 1 


1 °! 


CT\ 


<j\ 


<J\ 


G\ 


a\ 




cri 


<j\ 




<J\ 





a 3 o 



3 



fM 



fll 



QJ 



CN 






<D 




<3* 


r- 


<J\ 


O 


CN 


ID 


r- 


in 


in 


<j\ 


CN 






00 


o 


ro 


00 


rH 


in 


in 


O 


ro 


VD 


rH 


CN 


^D 






CN 


in 


ro 


rH 




in 




rH 


a\ 


CN 


-^r 


in 


CN 


1 


1 


CN 


CN 


o 


rH 


rH 




o 


rH 




O 


o 


ro 




1 


1 






rH 




















rH 






r^ 


in 


CT» 


rH 


O 


ro 


ID 


<J\ 


ro 


00 


o 




00 








ro 


r- 


00 


O 


ro 


00 


ro 


o 






CN 


o 






o 


i — 1 


ID 


00 


in 


00 


in 


ID 


00 


CN 


'vD 


00 


• 


1 




• 


• 


• 


• 


» 


• 


• 


• 


,o 


• 


• 


• 


in 


1 


1 






O 


in 




rH 




ro 


t 






o 


rH 






00 


<T> 


<J\ 


in 


cr» 


ro 


o 


VD 


00 


VD 


rH 


r- 


ro 






ID 


ro 


in 


ro 


CN 


ro 


a\ 


in 


CN 




00 






1 


1 




in 




ro 






o 


rH 


rH 


O 


o 


o 


o 


in 


in 


o 


in 


in 


o 


in 


in 


o 


in 


in 


O 


in 


in 


in 


CN 


o 


in 


CN 


o 


in 


CN 


o 


in 


CN 


o 


in 


CN 


o 


























in 


in 


in 


00 


00 


00 


in 


in 


in 


ro 


ro 


ro 


i — l 


i — 1 


rH 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


in 


in 


in 


in 


in 


in 


in 


in 


in 


in 


in 


in 


in 


in 


in 


rH 


rH 


rH 


rH 


rH 


i — l 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


i — l 


rH 



W 



g 



cn ro 



in 



vd 



oo 



<T> 



CN 



ro 



in 



20 



°. i 

0 H 






LD 



m 



m 



CN 



CN 



cn 



m 



CN 



CN 



cn 



0. 



CD 


O 


o 


cn 


CO 


<Tl 


00 


ID 


VO 




VO 


cn 


LD 


e'- 


n 


c m 


o 


o 


cn 


<n 


<n 


CT» 


cn 


<n 


cn 


00 


00 


<n 


en 


r* 


Cm S 


* 


• 


<n 


<n 


cn 


<T» 


cn 


<T» 


<T» 


Cn 


cn 


<n 


<n 


<n 




r — 1 




























aj ^0 


VO 




vo 


CN 


in 


P0 


o 




VO 


CN 


<n 










rH 




no 


r- 


O 


CN 




cn 


<n 


00 


i — I 


o 


CN 


q 3 


VO 


i— 1 


o 


m 


o 


CN 


00 


CN 




• <n 


<n 


CN 


in 


cn 


-H 


CN 


cn 


CN 


m 


• 


H 


CN 


00 


o 


O 


l n 


O 


o 


n 


Cm 


• 


• 


t 


• 


iH 






















iH 


rH 


CO 


in 


cn 


in 


CN 


in 






o 


H 




H 




LD 


H 


cn 


cn 


in 


m 


00 


00 


cn 


cn 


r- 


00 


H 


00 




r — 1 


r* 


i ^ 


cn 


VO 


00 


cn 


r* 






LD 


CN 






■K 


• 


• 


• 


• 


• 


* 


• 


• 


• 


9 


• 


• 


• 


• 




in 




1 ^ 


CN 


o 


in 


CN 


o 


LD 


CN 


o 


VO 


cn 


o 




CN 


o 


CO 


<n 


in 


r* 


i — I 


cn 


cn 


VO 


CN 


<T» 


00 


00 




in 




rH 


CN 


00 


cn 


H 




o 


o 


cn 


in 


VO 




■K 




C 30 


in 


in 




cn 


cn 


cn 


H 


1 — 1 


o 


o 


o 


o 



3 

•H M-t 
H W 



o 


o 


r* 


e'- 


00 


00 




vo 




vo 


cn 


n 


vo 


n 


o 


o 




en 


cn 


cn 


cn 


cn 


cn 


00 


00 


cn 


CO 


r* 


• 


• 


CTr 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


. cn 


cn 


cn 



(0 





vo 


cn 






cn 


cn 


vo 


cn 


cn 


vo 


cn 


H 


00 


LD 




0 


i ^ 


CN 




00 


vo 


o 


r* 


H 


00 




vo 


00 


H 


VO 


CN 


0 




<n 




cn 


00 


cn 


H 


O 


o 




cn 


vo 


CN 




cn 




CN 

• 


CN 

» 


1 


CN 


CN 


cn 


H 


cn 


00 


o 


o 


LD 


O 


O 


in 






r* 




VO 




00 


H 


00 


vo 




vo 


CN 


cn 


H 


vo 




00 


o 






VO 


00 


rr 


r- 


H 


00 




CN 


vo 


CN 


00 




n 


H 




r* 


00 


00 


H 


o 


00 


cn 


<n 


vo 


o 








• 


• 


| 


• 


» 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 




n 


n 


1 


cn 


cn 


o 


VO 


CN 


O 




CN 


o 


VO 




o 




r* 


o 




cn 


CN 


cn 


cn 


cn 


CN 


vo 


cn 


cn 


00 


vo 


cn 




. vo 


n 






cn 


in 




00 


in 


cn 


H 


cn 


in 


00 






r* 


CO 


1 




in 


in 


cn 


CN 


cn 


o 


H 


o 


o 


o 


o 


■&| 


j • 


• 


1 
























• 




o 


n 


n 


o 


in 


in 


o 


in 


in 


o 


in 


in 


o 


in 


LO 




n 


CN 


o 


in 


CN 


o 


in 


CN 


o 


in 


CN 


o 


in 


CN 


O 




























in 


in 


in 




00 


CO 


00 


in 


in 


in 


m 


cn 


cn 


H 


H 


H 


o 


o 


o 




o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


O 


O 


o 


o 


o 




o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


c 


cn 


cn 


cn 


cn 


m 


m 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


• 

Q 


<0 


I s * 


CO 


cn 


o 


H 


CN 


cn 




in 


vo 


r- 


00 


cn 


o 




H 


H 


H 


H 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


cn 



21 



. JJ 



cd 



cn 



cn 



cn 



cn 



CN 



CN 



cn 



CN 



CN 



cn 



cn 



CN 



03 


<n 


O 


0 


cn 


00 


cn 


vo 


CN 


vo 


in 


nr 


VO 


rH 


00 


nr 


4 H 


<T\ 


O 


0 


cn 


C n 


n 


<n 


<n 


<n 


<n 


00 


00 


cn 


p* 


p- 


•h m 


<n 


O 


0 


cn 


<T\ 


n 


cn 


<T\ 


<T\ 


<T\ 


cn 


cn 


cn 


a> 


cn 


(h w 




rH 


rH 


























rH 


cn 


nT 


m 


00 


rH 


p- 


vo 


rH 


00 


in 


p- 


00 


cn 


CN 


00 




p* 


VO 


cn 


C T\ 


LD 


rH 


in 


00 


00 


LD 


vo 


p- 


nr 


vo 


00 


§ -0 


<n 


rH 


00 


<n 


LD 


00 


00 


VO 


rH 


nr 


0 


00 


CN 


in 


rH 


•H 3 




ID 


p* 


rH 


00 


CN 


rH 


CN 


O 


0 


rH 


vo 


0 


0 


in 


(h 


• 




rH 






rH 






rH 
















0 


VO 


CN 


nr 


VO 


nr 


00 


p- 


cn 


rH 


cn 


VO 


rH 


CN 


p- 




0 




CN 


00 


nr 


cn 


cn 


vo 


00 


00 


CN 


rH 


CN 


p- 


00 




cn 


cn 


00 


cn 


vo 


LD 


rH 


rH 


m 


vo 


nr 


LO 


vo 


nr 


nr 


■K 

u 


CN 


CN 


d 


nr 


CN 


d 


nr 


CN 


d 


LO 


CN 


d 


nr 


CN 


0 




O 


rH 


00 


00 


00 


vo 


rH 


<T\ 


0 


00 


00 


0 


00 


O 


vo 




, CO 


nr 


vo 


0 


cn 


00 


vo 


nr 


cn 


rH 


cn 


0 


nr 


in 


nr 


* I 


VO 


p* 


vo 


m 


nr 


m 


CN 


CN 


CN 


rH 


0 


rH 


0 


0 


0 


a} 

































■a 

fa 



00 


cn 


0 


cn 


00 


00 


vo 


CN 


VO 


in 


cn 


LD 


rH 


00 


nr 


cn 


cn 


0 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


p* 


00 


cn 


p- 


p* 


cn 


cn 


0 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 



rH 



‘3 

H- 





LD 


p* 


vo 


LD 


ID 


CN 


nr 


nr 


P- 


cn 


rH 


vo 


nr 


00 


00 


,0 


VO 


VO 


cn 


CN 


P* 


CD 


00 


0 


O 


ID 


P- 


cn 


nr 


cn 


nr 


3 


vo 


00 


p* 


rH 


CN 


O 


CD 


p- 


CD 


nr 


CN 


ID 


CN 


LD 


cn 




m 


LD 


LD 


CN 


CD 


O 


rH 


CN 


cn 


O 


rH 


vo 


0 


O 


nr 








rH 






rH 






















CO 


vo 


VO 


rH 


CN 


O 


' O 


CN 


vo 


P* 


nr 


ID 


p- 


cn 


00 




m 


rH 


cn 


VO 


ID 


O 


O 


ID 


ID 


00 


cn 


ID 


cn 


nr 


rH 




00 


cn 


nr 


nr 


O 


cn 


O 


O 


vo 


nr 


vo 


LD 


cn 


rH 


LD 


1 U 


rH 


rH 


d 


nr 


CD 


d 


nr 


CN 


d 


ID 


rH 


O 


^r 


CN 


O 




CN 


O 


rH 


rH 


cn 


CN 


CN 


cn 


CN 


0 


O 


VO 


rH 


CD 


cn 




cn 


O 


CD 


CO 


CN 


00 


ID 


CD 


CD 


rH 


P- 


O 


ID 


nr 


nr 




vo 


P* 


P' 




LD 


vo 


CN 


CN 


CD 


rH 


O 


rH 


O 


O 


0 


' cu 


• 


• 






























0 


LD 


LD 


O 


LD 


ID 


O 


ID 


ID 


O 


ID 


ID 


O 


LD 


LD 


5H 


IT) 


CN 


O 


LD 


CN 


O 


ID 


CN 


O 


ID 


CN 


O 


ID 


CN 


O 




























ID 


LD 


LD 




00 


00 


00 


LD 


ID 


ID 


CD 


CD 


CD 


rH 


rH 


rH 


O 


O 


O 


a 


0 


O 


O 


O 


O 


d 


d 


d 


d 


d 


d 


O 


d 


d 


O 


c 


0 


O 


O 


O 


O 


0 


0 


O 


0 


0 


0 


O 


0 


O 


O 




LD 


LD 


LD 


LD 


LD 


ID 


ID 


ID 


ID 


ID 


ID 


ID 


ID 


LD 


LD 




rH 


CN 


CD 


nr 


ID 


VO 


P- 


00 


cn 


O 


rH 


CN 


CD 


nT 


LD 


2 


cn 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


nr 


nr 


nr 


nr 


nr 


nr 



22 



4-1 

O 




co CN CN 



co cn 



CO CN CN 



CM CM CM 



CN CN CM 



1 «« 
•H <P| 

W 



o o o 



o 


o 


o 


G\ 




cn 


r^ 




r- 




r- 


LD 


CO 




rH 


o 


o 


o 






cn 


a>i 


a>i 


a>> 


a>i 


00 


00 


<n 


00 




• 


• 


• 




cn 


cn 


<j\ 


a>i 




a>i 


a\ 


cn 


cn 


CTi 








1— 1 



























i— i 


00 


O 


CN 




o 


CN 


00 


CN 


rH 


o\' 




rH 




VO 


VO 


Crf 


in 


<sT 


O 


VO 


CN 




<n 


in 




rH 


o 


r^ 


CO 


CN 


CO 


c o 


VO 




CN 


00 


CN 


VO 


CM 




CO 


<nT 


o 


00 


CN 


LO 


CO 


■H 3 


! CM 






rH 


CO 


CN 


rH 


CN 


CN 


o 


rH 


VO 


O 


o 




(■H 






rH 






rH 






rH 
















o 


rH 


CN 


rH 


G\ 


LD 


o 




00 


00 


o 


VO 


VO 


VO 


<n 




in 


CO 


00 


o 




CN 


r- 




00 


CN 


CN 


a\ 




e'- 






in 


00 


in 




CN 


in 


VO 


LD 


CO 


in 


00 




in 


en 


LD 


* 


• 


• 


• 


■ 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


P 


in 


CN 


o 


in 


CO 


o 




CN 


o 


in 


CN 


o 


LD 


CN 


O 




in 


00 


VO 


CO 


in 




CO 


CO 


CO 


rH 


CN 


rH 


VO 


<n 


O 




CM 


CN 


VO 


rH 


VO 


o 


00 


00 


VO 


O 


rH 


<n 


LD 


LO 


in 


& 


00 


00 


00 


in 


in 


in 


CN 


CN 


CN 


rH 


rH 


o 


O 


O 


o 


rH 
































nJ 
































■H 


o 


O 


O 


























-P 4H 


o 


O 


O 


00 


00 


<n 


VO 








[■"> 




CO 


<vT 


o 


■a a 


o 


O 


O 


G\ 


<n 


G\ 


G\ 


o\ 


cn 


cn 


00 


00 


cn 


00 


r- 


• 


• 


• 




o\ 


<n 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


H 


rH 


rH 


rH 



























•rH 

4-> 




? P 



? a 



u 

CU 

c 








VO 


o 


CO 


O 


CM 


<n 


00 


r- 


CO 


o 


CN 


CO 


CO 


r- 


00 


CO 


r- 


00 




O 


' <vT 


in 


<n 


CM 


o 




CO 


CN 


<vT 


in 


VO 


CO 


O 


in 


00 


CN 




CO 




o 




CN 


in 


cn 


CN 


<vT 


• 


CN 


CO 


CM 


rH 


CN 


CN 


o 


rH 


VO 


o 


o 


CO 


« 


• 


rH 



























00 ' 


in 


VO 


VO 


00 


o 


in 


CO 




VO 


r- 


CO 


CO 


o 


r- 


in 


CN 


CN 


rH 


00 


rH 


rH 


CO 


00 


in 




in 


CN 


CN 


in 


r- 


<n 


VO 


VO 


["- 


in 


rH 


in 


CO 


CO 


00 


in 


VO 


o 


vo 


in 


CN 


d 




CN 


o 


in 


CN 


o 


in 


CN 


o 


in 


CO 


o 



o 


CN 


■^r 


CO 


in 


o 


CN 


CN 


rH 


00 


CO 


O 


VO 


o 


r- 


CO 


CO 


r- 


vo 


CN 


o 


o 


00 


vo 


cn 


rH 


O 


in 


vo 


in 


00 


00 


00 




LO 


in 


CO 


CN 


CN 


o 


rH 


rH 


o 


o 


o 



o 


in 


in 


o 


in 


in 


o 


in 


in 


o 


in 


in 


o 


in 


in 


in 


CN 


o 


in 


CN 


o 


in 


CN 


o 


in 


CN 


o 


in 


CN 


o 


























in 


LO 


in 


00 


00 


00 


in 


in 


in 


CO 


CO 


CO 


rH 


rH 


rH 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 



o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


VO 




CO 


<n 


o 


rH 


CN 


CO 




in 


vo 


r- 


CO 


CTi 


o 










in 


in 


in 


in 


in 


in 


LO 


in 


LO 


in 


vo 



23 



IV. RESULT AND CONCLUSION 



In this chapter efficiencies for various p-r pairs were 
given in Table (4.1). The method developed was used to 
calculate table entries. Three entries were given for each 
p-r pair. Entry in the middle is efficiency value for that 
pair at optimum t or u Q . Entries inside parenthesis are 
also efficiencies and found as follows: 

Let u = - log t 
o o 

If we miss u Q by ±50 percent, what would the efficiency be? 
Efficiency inside the parenthesis above the optimum efficiency 
answers the question for u-^ = 0.5 u Q , and other one below 
the optimum answers for U£ = 1.5 u q . 

The change u = -log t (or rather its inverse t = e U ) 

X 

changes the probability generating t function into the 
empirical La Place transform. It is seen that a 50 percent 
error in its argument u results in only a minor degradation 
of efficiency, whereas a 50 percent error in t would have 
disasterous effects. In all of the examples, the initial 
u Q was well within this 50 percent range, the worse case 
being 38 percent. 

After examining Table (4.1) one can conclude that effi- 
ciency is monotonically increasing with increasing p and 
increasing r for r greater than one and minimum efficiency 
at optimum t or u Q is greater than 0.95. 
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Following Table (4.1) optimum t and u Q (in parenthesis) 
were tabulated in Table (4.2). One may observe that t 

o 

decreases with p and increases with r, whereas u q behaves 
oppositely. 

After Table (4.2) efficiency contours were graphed in 
p-r plane (Figure 4.1), in mean-r plane (Figure 4.2) and 
in mean-exp (r) plane (Figure 4.3). Compared with the 
Figure in [5] Figure (4.3) always gives better efficiency 
values than all the other alternatives. 
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TABLE 4.1. ASYMPTOTIC EFFICIENCY OF p* , r* (Cont'd) 
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TABLE 4.1. ASYMPTOTIC EFFICIENCY OF p* , r* (Cont'd) 
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TABLE 4,1, ASYMPTOTIC EFFICIENCY OF p * , r* (Cont'd) 
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981) (.981) (.981) (.982) (.982) (.982) (.984) 



TABLE 4,1. ASYMPTOTIC EFFICIENCY OF p* , r* (Cont’d) 
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TABLE 4.1. ASYMPTOTIC EFFICIENCY OF p* , r* (Cont'd) 
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TABLE 4.1. ASYMPTOTIC EFFICIENCY OF p* , r* (Cont'd) 
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TABLE 4,1, ASYMPTOTIC EFFICIENCY OF p* , r* (Cont'd) 
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TABLE 4,1. ASYMPTOTIC EFFICIENCY OF p*, x* (Cpnt'd) 
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TABLE 4.1. ASYMPTOTIC EFFICIENCY OF p* , r* (Cont'd) 
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Figure 4.3 Efficiency contours in MEAN-EXP (r) plane 
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APPENDIX A 



DETAILED COMPUTATIONS 

The two estimating equations g^ and g 2 given in chapter 
II are 



STi = 



x - sa = o 



- / p x r 

q 2 t ~ 1-qt 



= 0 



Matrix A's entries were given as follows: 



11 



39 i 

- E( 3r> * 



E(^) 



12 



8g. 

= E (^) = 

3r 



E(- §) 



q 

p 



21 



3g 9 

B ^ 1 



1 qt (1-qt) 



- -<r5t> 



p v r-1 r (1-t) 



(1-qt) 



22 



ag 9 

E <inr> 
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Matrix C's entries were given as follows: 



C 



11 



n E(g 2 ) = n E[(x - ^) 2 ] 

n 

n Var(x) = n Var(^- £ x . ) 

n u 1 

i=l 



n 



n— ir Var ( \ x. ) 
n . , 1 



i=l 



£a 

2 

p 



Var (x) 



C 12 = n E lV 9 2 ] 

_ — , n — 

= n Cov (X, t x ) = n Cov (— \ x.,t x ) 

i=l 

= Cov ( x if t x ) 

n n n 

= l Cov(x i ,t X ) = l Cov(x i , 

i=l i=l j=l 

n n 

= ~ l l Cov(x i ,t Xj ) 

i=l j=l 
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»H | G 



% 



X . 

Since Cov(x^,t =0 for i ^ j because of independence, 
it follows that 



n 



C 12 C 21 



— y cov (x . , t 1 ) 
n L 1 

i=l 



= Cov (x, t x ) = E(x-t x ) - E(x)E(t x ) 



One can calculate E(x*t ) using the following relation. 



G ( t) = E(t X ) G' (t) = E(xt X-1 ) 



Then 



E(xt X ) = tG' (t) 



where 



G(t) 



i - qt 



r 



taking the derivative and multiplying by t, E(xt x ) is found 
as 



E (xt X ) 



trqp r 

(1-qt) 



r+1 
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Then it yields 



r _ r . rqp r 1 (l-t) 

12 " 21 "(l-qt) r+1 

C 22 > n E(g 2 ) = El(t x - (j^) 1 ) 2 ] 

n 

.1 . . ^ 

= n Var (t X ) = n Var(i Ft 1 ) 

n L 

i=l 

= £ Var( l t x ) 

i=l 

= Var (t x ) = E (t 2x ) - [ E { t x > ] 2 



G(t 2 ) - [G(t)l 2 



<-E-j) r - <T^7F) 2r 
1-qt 2 1 " 3t 



As the reader will recall, | A | was given by (2.4) and can be 
written as follows. 



m 



- -V i 



n 



r!n! 



n 



n+1 (n+r) ! 



_ r In I 
n+1 (n+r) ! 



•] 



n=l 



n=m+l 



Call the summation from n = m+1 to «> the remainder . The 
remainder can be adjusted such that it will be smaller than 
any specified small value (e) . 
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n 



r ! n! 



n 



n+1 (n+r) ! — 



n+1 



n=m+l 



n=m+l 



r!n! 



Because ( n + r ) i £ 1 (reciprocal of binomial coefficient) , 



.n 



m+1 

3 J g» 



g 



m+1 



^ n+1 — m+2 

n=m+l n=0 



(m+2)p - 



< e 



A value of m can be found such that 



m+1 



(m+2) p - 



< e 



holds. Then this m guarantees the remainder will be smaller 
than e, and the efficiency equation of (2.12) becomes 



Eff = £E- 



r+2 



r l-qt , , 1-qt . _ q(l-t) . 2 

P P P 



t ( 1-qt) 



2 r 
(1+qt: ) r 



2(r+1) r„ .,2 r„ ..2, ,?q“ ilnl , 
p (1-qt) -rqp (1-t) ] ( 2, r-rr rzrrzn . ) 



1 '•.n+l (r+n) ! 
n=l 
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APPENDIX B 



SAMPLES GENERATED BY COMPUTER 

Samples generated by computer for sixty cases used in 
chapter three, example four, are given in the following 
pages case by case. 

See subroutine GNRNB of computer program- 3 of Appendix 
C for the simulation procedure, which was taken from [6] . 
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O i— I CN ro LO O H CM o OHCMn^lOhC^OJ o H cm n M* ui r- O i — I CM to 



VAL 



VAL. 



FRQ . 



FRQ. VAL. FRQ. VAL. FRQ. 



CASE 1 



2 

7 

1 

3 

2 

CASE 2 



8 

5 

2 

CASE 3 



15 

CASE 4 



1 

1 

1 

4 

1 

3 

2 

1 

1 

CASE 5 



2 

1 

4 

1 

2 

1 

4 

CASE 6 



8 

5 

1 

1 





CASE 


7 




CASE 


10 




CASE 


13 










(Cont'd) 




(Cont ' d) 


3 




2 














5 




1 


63 




1 


68 




1 


8 




1 


64 




1 


70 




1 


9 




1 


71 




1 


88 




1 


10 




2 


73 




1 


120 




1 


11 




1 


74 




1 


131 




1 


12 




1 


76 




1 


139 




1 


13 




1 








153 




1 


14 




1 




CASE 


11 


176 




1 


15 




1 














19 




1 


2 




1 




CASE 


14 


21 

22 




1 

1 


7 

8 




1 • 
1 


12 




1 








9 




1 


15 




1 




CASE 


8 


11 




1 


24 




1 








17 




1 


30 




1 


1 




1 


18 




1 


40 




1 


2 




2 


19 




1 


42 




1 


3 




3 


22 




2 


50 




1 


5 




1 


25 




1 


52 




1 


6 




1 


27 




1 


59 




1 


7 




1 


30 




1 


62 




1 


8 




3 


37 




1 


70 




1 


9 




1 


41 




1 


71 




1 


11 




1 








78 




1 


14 




1 




CASE 


12 


85 




1 














96 




1 




CASE 


9 






A 














u 

1 




2 




CASE 


15 


0 




5 


4 




1 








1 




2 


5 




1 


0 




3 


2 




2 


6 




3 


1 




2 


4 




2 


7 




1 


2 




1 


5 




2 


10 




1 


3 




1 


6 




2 


11 




1 


5 




1 








25 




1 


9 




2 




CASE 


10 








15 




1 










CASE 


13 


18 




1 














22 




1 


11 




1 








30 




1 


14 




1 


33 




1 






l 


21 




1 


49 




1 


3 ** 




X 


30 




1 


52 




1 




CASE 


16 


33 




2 


59 




1 








39 




1 


60 




1 








40 




1 


63 




1 


0 




8 


57 




1 


64 




1 


l 




7 



2 6 

3 6 

4 1 

5 2 
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VAL 



FRQ. 

CASE 17 



VAL. FRQ. 
CASE 22 



VAL. FRQ. 
CASE 25 



VAL. FRQ. 
CASE 27 



0 14 

1 8 

2 5 

3 3 

CASE 18 



0 27 

1 3 

• CASE 19 



1 2 

2 6 

3 7 

4 3 

5 4 

6 3 

7 1 

10 1 

11 2 

14 1 



CASE 20 



0 

1 

2 

3 

4 

5 
7 



5 

8 

2 

2 

7 

4 

2 



CASE 21 



0 

1 

2 

3 

4 



19 

5 

3 

2 

1 



1 1 

3 2 

6 1 

7 3 

8 4 

9 2 

11 5 

12 1 

13 1 

14 2 

15 1 

16 1 

17 1 

19 2 

20 2 

26 1 



CASE 23 



0 3 

1 2 

2 2 

3 3 

4 4 

5 5 

6 4 

7 1 

8 2 

10 2 

15 1 

20 1 



CASE 24 



0 

1 

2 

3 

4 

5 
7 



12 

10 

2 

1 

1 

2 

2 



14 

23 

28 

29 

30 

31 

32 

33 

34 

36 

37 
47 
51 
53 
58 
60 

67 

68 
71 
73 
77 
81 
88 

102 



1 

1 

3 
1 
1 
1 
1 

4 
1 
1 
2 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



CASE 26 



0 1 

5 1 

6 1 

8 1 

9 1 

10 1 

12 1 

15 3 

16 2 

17 1 

20 2 

21 1 

23 1 

24 1 

27 2 

29 1 

31 1 

32 1 

33 1 

35 1 

47 1 

49 1 

53 1 

61 1 



0 8 

1 4 

2 1 

3 3 

4 3 

5 2 

7 1 

8 1 

9 2 

11 1 

15 1 

18 1 

26 1 

29 1 



CASE 28 



42 1 

46 1 

52 1 

55 1 

57 1 

58 1 

59 1 

60 1 

65 1 

72 1 

73 1- 

79 1 

81 1 

91 1 

92 1 

98 1 

100 1 

106 1 

112 1 

118 1 

120 1 

121 1 

124 1 

129 1 

132 1 

141 1 

148 1 

167 1 

185 1 

197 1 
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VAL 



FRQ 



FRQ. 



VAL . FRQ . 



VAL . FRQ . 



VAL. 



CASE 29 



CASE 31 



CASE 36 



CASE 39 



7 


1 


0 




21 


12 


1 


1 




13 


15 


1 


2 




11 


18 


1 


3 




4 


21 


3 


7 




1 


22 


1 








26 

30 


1 

1 




CASE 


32 








33 

40 

49 

52 

54 

55 
57 


1 

1 

1 

1 

4 

2 

1 


0 

1 

2 

3 

4 

5 




24 

17 

6 

1 

1 

1 


58 

65 


1 

1 




CASE 


33 


66 


1 








68 


■ 1 


0 




44 


71 


2 


1 




4 


74 


1 


2 




2 


86 


1 








95 


1 




CASE 


34 


CASE 


30 


1 




8 






2 


* 


6 


0 


6 


3 




3 


1 


5 


4 




7 


2 


2 


5 




10 


3 


2 


6 




4 


4 


1 


7 




2 


5 


2 


8 




4 


6 


1 


9 




4 


8 


1 


11 




1 


9 


1 


17 




1 


10 


1 








12 

15 


2 

1 




CASE 


35 








21 

25 


1 

1 


0 

1 




9 

9 


38 

47 


1 

1 


2 

3 




8 

8 


52 


1 


4 




5 






5 




6 






6 




1 






7 




2 






8 




1 






9 




1 



0 

1 

2 

3 



37 

6 

6 

1 



CASE 37 



2 1 

3 3 

4 3 

5 2 

6 2 

7 2 

8 5 

9 3 

10 3 

11 2 

12 6 

13 1 

14 2 

15 4 

17 4 

20 1 

22 2 

23 1 

25 1 

29 1 

33 1 



CASE 38 



0 1 

1 7 

2 4 

3 3 

4 5 

5 6 

6 4 

7 6 

8 2 

10 3 

11 2 

14 1 

15 2 

16 2 

17 1 

25 1 



0 

1 

2 

3 

4 

5 

6 

10 



27 

5 

8 

4 

.3 

1 

1 

.1 



CASE 40 



8 1 

16 1 

17 1 

22 1 

26 3 

27 4 

28 1 

30 1 

31 1 

32 3 

33 3 

35 2 

36 2 

37 1 

40 2 

41 1 

44 2 

45 1 

47 1 

48 1 

49 1 

54 2 

55 1 

58 1 

60 1 

61 1 

62 1 

68 1 

70 1 

71 2 

74 1 

75 1 

87 1 

91 1 

93 1 
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VAL. FRQ. 


VAL 


FRQ 




CASE 41 




CASE 43 


3 


1 


17 


1 


4 


2 


21 


1 


8 


3 


33 


1 


9 


3 


40 


1 


10 


3 


43 


1 


11 


1 


. 45 


1 


12 


1 


50 


1 


13 


1 


51 


2 


14 


4 


52 


1 


15 


4 


55 


1 


16 


1 


57 


1 


17 


2 


61 


1 


18 


3 


63 


2 


19 


1 


64 


1 


20 


2 


68 


1 


21 


1 


69 


1 


22 


1 


71 


1 


23 


1 


76 


1 


24 


1 


78 


1 


26 


1 


80 


1 


27 


3 


82 


1 


29 


1 


83 


1 


33 


1 


87 


1 


35 


1 


90 


1 


41 


1 


93 


1 


47 


1 


94 


1 


50 


2 


95 


1 


74 


1 


98 


1 


76 


1 


100 


1 


82 


1 


101 


1 






102 


1 




CASE 42 


104 


1 






109 


1 


0 


16 


114 


2 






121 


1 


X 

9 


Q 


125 


1 


z 


y 

9 


127 


1 




Z 

i 


129 


1 




X 


135 


1 


D 

£ 


9 


140 


1 


D 

n 


Z 


149 


1 


/ 

o 


X 

9 


158 


1 


o 

10 


2 


163 


1 


12 


3 


167 


1 


16 


2 


171 


1 


20 


1 


184 


1 


35 


1 


185 


1 
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APPENDIX C 



COMPUTER PROGRAMS 

Computer programs, to prepare Table (4.1) and Table (4.2) 
(Computer program 1) , to find the estimates for the examples 
one through three of chatper three (computer program 2) , 
and to generate the data for the sixty cases of the example 
four of chapter three and to find the estimates for them 
(computer program 3) , are given in the following pages. 

There are two subroutines, namely GOLDEN (Golden section 
search), SVS (Single variable search), and one function, 
namely FCN (to evaluate efficiency function) , which are 
common to the three programs mentioned above. These are 
given at the beginning under the title COMMON SUBROUTINES. 
After that the three programs are given. 

One who wants to use any one of the programs, must include 
the common subroutines after the main program of the computer 
program he chooses. 
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COMMON SUBROUTINES 



SUBROUTINE S V S ( FCN , TNOT ) 

IMPLICIT REAL*8 (A-H,0-Z) 

COMMON P ,R, DETLAM 
Q=l. ODO-P 

UPPER=1 . ODO/DABS (Q) 

ALOW=-l . ODO/DSQRT (DABS (Q) ) 

RIOU=0 . 000001DO 
START=0 .9999 
U= START 
S=0. 0001D0 
A=FCN(U) 

Ul=U+S 
B=FCN (Ul) 

IF (B.GT.A) GO TO 5 
IF (B.EQ.A) GO TO 6 
M=3 

1 U2=START+S*2 . 0D0** (M-2) 

IF (U2.GE. UPPER) GO TO 99 
A=B 

B=FCN (U2 ) 

IF(A.EQ.B) GO TO 7 

IF (A.LT.B) GO TO 8 

M=M+1 

U=U1 

U1=U2 

GO TO 1 

7 CALL GOLDEN (FCN, Ul , U2 , RIOU , 50/TNOT) 

GO TO 999 

99 CALL GOLDEN (FCN ,Ul /UPPER /RIOU/ 50/TNOT) 

GO TO 999 

8 CALL GOLDEN (FCN, U , U2 , RIOU , 50 , TNOT) 

GO TO 999 

6 M=2 

CALL GOLDEN (FCN , U , Ul , RIOU , 5 0 , TNOT ) 

GO TO 999 
5 M=3 

2 U2=START-S* 2 . 0D0** (M-3 ) 

IF (U2.LE.ALOW) GO TO 77 
B=A 

A=FCN (U2) 

IF (A.EQ.B) GO TO 17 

IF (A.GT.B) GO TO 18 

U1=U 

U=U2 

M=M+1 

GO TO 2 

17 CALL GOLDEN (FCN, U2,U, RIOU ,50, TNOT) 

GO TO 999 

18 CALL GOLDEN (FCN, U2,U1, RIOU, 50, TNOT) 

GO TO 999 

77 CALL GOLDEN (FCN, ALOW, U, RIOU, 50, TNOT) 
99? RETURN 

END 
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SUBROUTI NE GOLDEN ( FCN , ENDL , ENDR , RIOU , NMAX , XOLD ) 
IMPLICIT REAL* 8 (A-H,0-Z) 

COMMON P , R, DETLAM 
SIGMA= ( 3 . ODO-DSQRT (5 . ODO) )/2.0D0 
AL=ENDR-ENDL 
— - XOLD=ENDL+SIGMA*AL 
FOLD=FCN (XOLD) 

DO 1 1=2, NMAX 
XNEW=ENDL+ENDR-XOLD 
FNEW=FCN (XNEW) 

IF (FNEW.LE.FOLD) GO TO 5 

IF (XNEW. LT. XOLD) GO TO 6 

ENDR=XNEW 
20 AL=ENDR-ENDL 

GO TO 30 

6 ENDL=XNEW 
GO TO 20 

5 IF ( XNEW. LT. XOLD) GO TO 7 

ENDL=XOLD 
8 XOLD=XNEW 

FOLD=FNEW 
GO TO 20 

7 ENDR=XOLD 
GO TO 8 

30 IF (AL.LE.RIOU) GO TO 40 
1 CONTINUE 

40 RETURN 

END 



FUNCTION FCN (X) 

IMPLICIT REAL *8 (A-H,0-Z) 

COMMON P,R, DETLAM 
ZR = 1.0D-73 
Q=l. 0D0-P 
A=R*P** (R+2 . DO) / Q 
u-" B=1.0D0-Q*X 

C= (DABS ( (B*DLOG (B/P)/P-Q* (1.0D0-X)/P) ) ) **2.D0 
D=B** ( 2 . D0*R+2 . 0D0) / (1 . 0D0-Q*DABS (X) **2.D0) **R 
E=P**R*B**2.D0+R*Q*P**R* (DABS (1.0D0-X) ) **2.D0 
G=-A*C/DETLAM 
F=D-E 

IF (DABS (F) .LT.ZR) GO TO 5 

FCN=G/F 

GO TO 7 

5 FCN=-1 . ODO 

7 RETURN 

END 
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COMPUTER PROGRAM 1 



IMPLICIT REAL* 8 (A-H,OZ) 

EXTERNAL FCN 

DIMENSION PR (14) ,RR(26) ,EFF(26 / 14) ,TZERO(26,14) ,TIEST(52) ,EFF1(26, 
114) / EFF2 (26,14) 

COMMON P,R,DETLAM 

DATA EFF , TZERD/728*0 . 0D0/ 

DATA EFF1,EFF2/728*0.0D0/ 

DATA PR/.05D0, .06D0, .07D0, .08D0, .09D0, .IDO, .2D0, .3D0, .4D0, .5D0, .6D 
10,.7D0,.8D0, .9D0/ 

DATA Riy'.5DO,l.DO,1.5DO,2.DO,2.5DO,3.DO,3.5DO,4.DO,4.5DO,5.DO,6.DO 
1,7. DO, 8. DO, 9. DO, 10. DO, 11. DO, 12. DO, 13. DO, 14. DO, 15. DO, 16 .DO, 17 .DO, 18 
l.D0,19.D0,0.0D0,0.0D0/ 

DATA TTEST/52*0.0D0/ 

WRITE (6,444) 

444 FORMAT ('I') 

EPS=0. 0001D0 
DO 1 J=1 , 14 
P=PR(J) 

Q=1 . ODO-P 
K=1 

7 TEST=Q**(K+1)/(P*(DFLQAT(K)+2.0D0)) 

K=K+1 

GO TO 7 

8 M=K 

WRITE (6,21) M 

21 FORMAT (15X, 'UPPER BOUND=' ,18) 

DO 2 1=1,24 
R=RR(I) 

DETLAM=0 . ODO 
DO 10 N=1,M 
RATI0N=1 . ODO 
DO 12 JJ=1,N 

12 RATICN=RATION*DFLQAT ( JJ) / (DFLQAT ( JJ) +R) 

DETLAM=DETLAM+0**N*RATI(^/DFLCAT (1SH-1) 

10 CONTINUE 

CALL SVS (FCN,TNOT) 

TZERO (I , J) =TNOT 

TPL=DEXP ( 0 . 5DO*DLOG (TNOT) ) 

mEN=DEXP (1 . 5DO*DLOG (TNOT) ) 

EFF2 (I , J) =-FCN (TPL) 

EFF1 ( I , J) =-FCN (TMIN) 

EFF (I, J) =-FCN (TNOT) 

2 CONTINUE 

1 CONTINUE 

WRITE (6,200) (PR(I) ,1=1,14) 

200 FORMAT (///,50X,'P VALUES' ,/,12X,14 (2X,F5. 3, IX) ,/,12X,14 (2X, ' 

1\1X),/) 

DO 13 1=1,24 

WRITE (6,202) (EFF1 (I,J) , J=l,14) 

WRITE (6,201) RR(I) , (EFF(I,J) ,J=1,14) 

WRITE (6,202) (EFF2 (I, J) , J=l,14) 
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201 



FORMAT (3X, 'R=' ,F5.2,2X,14 (2X,F5.3,1X) ) 
WRITE (6,203) (TZERO(I,J) ,J=1,14) 

202 FORMAT (12X, 14 (IX, ' (',F5. 3,') ')) 

203 FORMAT (5X, 'OP.TNOT' ,14 (2X,F5.3,1X) ,//) 

13 CONTINUE 

STOP 

END 
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COMPUTER PROGRAM 2 



IMPLICIT REAL* 8 (A-H,0-Z) 

EXTERNAL FCN 

DIMENSION FRQ(20) ,OBS(20) 

COMMON P , R, DETLAM 
WRITE (6,444) 

444 FORMAT ( ' 1 ' ) 

EPS=0 . 0001D0 
WRITE (6,91) 

91 FORMAT (5X, 'READ IN NO. OF OBS. AND VARIATES IN 215 FORMAT',/) 

READ (5,81) NOBS , . NVAR 

81 FORMAT (215) 

WRITE (6,92) 

92 FORMAT (5X, 'NOW GIVE ME VARIATE VALUES AND CORRESPONDING FREQUENCY 

UN THE FORMAT 2F10.5 FOR EACH CARD' ,/) 

SIMO.ODO 

SUM1=0 . 0D0 

DO 4 1=1, NVAR 

READ (5,82) CBS (I) ,FRQ (I) 

82 FORMAT (2F10. 5) 

SUMNSUM+OBS (I) *FRQ(I) 

SUM1=SUM1+ (OBS (I) **2. 0D0) *FRQ (I) 

XBAR=SUM/DFLQAT (NOBS) 

SSQR= (SUMl-DFLQAT (NOBS) *XBAR**2 . 0D0) /DFLOAT (NOBS-1) 

WRITE (6,209) XBAR, SSQR 

209 FORMAT (15X, 'USING GIVEN DATA FOLLOWING INFORMATION IS OBTAINED',// 
1,2 OX, 'SAMPLE MEAN =' ,F7. 4,//, 20X, 'SAMPLE VARIANCE =',F7.4,// 

2 ) 

P=XBAR/ SSQR 

R=XBAR**2 . 0D0/ (SSQR-XBAR) 

. Q=1 . ODO-P 
K=1 

7 TEST=Q** (K+l) / (DFLOAT (K+2) ) 

IF (TEST.LE.EPS) GO TO 8 
K=K+1 

GO TO 7 

8 M=K 
DETLAM=0 . 0D0 
DO 10 N=1,M 
RATIONAL. 0D0 
DO 12 JJ=1,N 

12 RATIQN=RATICN*DFIQAT(JJ)/ (DFLOAT (JJ)+R) 

DETIAM=DETIAM4^**N*RATION/D FLOAT (N+l ) 

10 CONTINUE 

1^1 

95 CALL SVS (FCN,TNOT) 

EFF1=-FCN (TNCT) 

RSTAR=R 
PSTAR=P 
SUM2=0 . 0D0 
DO 14 J=1,NVAR 

14 SUM2=SUM2+FRQ(J)*TNOT**OBS(J) 

TXBAR=SUM2/DFLQAT (NOBS) 
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CALL CQMPR(TNOT / TXBAR,XBAR,RSTAR) 

P=R/ (R+-XBAR) 

0=1 . ODO-P 
K=1 

77 TEST=Q* * (K+l) / (P^DFLQAl (K+2) ) 

IF (TEST.LE.EPS) GO TO 88 
K=K+1 
GO TO 77 
88 M=K 

DEILAM=0 . ODO , 

DO 51 N=1,M 
RATIONAL. ODO 
DO 52 JJ=1,N 

52 RATTON=RATION*DFLQAT ( JJ) / (DFLQAT ( JJ) +R) 

51 DETLAM=DETLAI^**N*RATION/DFIOAT (N+l) 

EFF2=-FCN (TNOT) 

IF (DABS (P-PSTAR) .GT. 0.005) GO TO 29 
IF (DABS (R-RSTAR) .GT. 0.005) GO TO 29 
IF (DABS (EFF1-EFF2) .GT. 0.04) GO TO 29 
WRITE (6,94) P,R,EFF2,TNCT,L 

FORMAT (20X, ’ESTIMATE OF P =’ ,F7.5,//, 20X, ’ESTIMATE OF R =’, 

17 . 4 , // , 2 OX , EFFICIENCY =’ ,F7. 5,//, 20X, ’OPTIMUM TNOT =’, 

27.5, //,20X, ’NO. OF ITERATIONS =’,FI4,//) 

GO TO 999 

29 L=L+1 

GO TO 95 
999 STOP 

END 



SUBROUTINE COMPR(T,TX,XB,RST) 

IMPLICIT REAL*8 (A-H,0-Z) 

COMMON P,R,DETLAM 
R1=RST 

30 B+R1+XB* (1 . ODO-T) 

TY= (Rl/B) **R1 
F=TX-TY 
FP1=-TY 

FP2= (XB* (1 . ODO-T) /B) +DL0G (Rl/B) 

FP=FP1*FP2 

R2=R1-F/FP 

IF (DABS (R2-R1) .LE.0. 00001D0) GO TO 66 
R1=R2 
GO TO 30 
66 R=R2 

RETURN 

END 
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COMPUTER PROGRAM 3 



IMPLICIT REAL* 8 (A-H,0-Z) 

EXTERNAL FCN 

DIMENSION FRQ(200) ,OBS(200) ,PVECT(10) ,RVECT(10) ,NVECT(10) ,DAT(200) 
1,ISD(200) 

COMMON P,R,DETLAM 

DATA RVECT/5 . 0D0 , 2 . 5D0 , 0 . 5D0 , 7*0 . 0D0/ 

DATA PVECT/ . 8D0 , . 5D0 , . 3D0 , . IDO , . 05D0 , 5* 0 . 0D0/ 

DATA NVECT/15, 30,50, 100 , 6*0/ 

CALL OVFLOW 
WRITE (5, 541) 

541 . FORMAT (//,15X, 'READ ISEED USED IN GENERATING DATA') 

READ (5,978) IS 

978 FORMAT (110) 

WRITE (6, 979) IS 

979 FORMAT (15X, 'ISEED USED TO GENERATE DATAf',110) 

CALL INT(IS,ISD,200) 

WRITE (6,444) 

444 FORMAT ( ' 1 ' ) 

EPS=0 . 0001D0 
NCASE=0 
DO 112 IN=1,4 
DO 113 IP=1,5 
DO 114 IR=1, 3 
P=PVECT(IP) 

R=RVECT(IR) 

NOBS=NVECT (IN) 

NCASE=NCASE+1 
NCP=2*NCASE 
NCM=NCP-1 
ISEED1=ISD (NCP) 

ISEED2=ISD (NCM) 

CALL GNRNB (NOBS , OBS , P , R, ISEED2 , ISDEED1 ) 

OBMX=-99999.DO 
DO 321 LK=1 , NOBS 

IF (QBS(LK) .GT.OBMX) OBMX=OBS(LK) 

321 CONTINUE 

OBIM=OBMX+10.DO 

NVAR=0 

390 OBMt^K>BMX+5.DO 

DO 370 MI=1 , NOBS 

IF (OBS (MI) .LT.OBMN) OBMN=0BS(MI) 

370 CONTINUE 
ICNT=0 

DO 380 MJ=l,NOBS 
IF (OBS(MJ) .NE.OBMN) GO TO 380 
ICNT=ICNTM-1 
OBS (MJ) =OBIM 
380 CONTINUE 

NVAR=NVAR+1 

DAT (NVAR) =OBMN 

FRQ (NVAR) =DFLQAT (ICNT) 

IF (OBMN.EQ.OBMX) GO TO 395 
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GO TO 390 

395 WRITE (6,889) 

889 FORMAT (//,10X, 89 ('*') ,/,40X, 'DATA GENERATED IS GIVEN BELOW' ,/,39X, 
131 ( '*' ) ,/) 

WRITE (6,301) (DAT(LO) ,L0=1,NVAR) 

301 FORMAT (20X, 'VALUES’ ,/,20X, 6 ('-') ,/, (5X,20 (1X,F5.0) ,/) ) 

WRITE (6,302) (FRQ(LO) ,L0=1,NVAR) 

302 format920X, 'FREQUENCIES ' ,/, 20X,11 ( '- ' ) ,/, (5X,20 (1X,F5. 0) ,/) ) 
WRITE (6, 543) 

543 FORMAT (//,10X, 'DATA WRITTEN ABOVE IS BASED ON FOLLOWING PARAMETERS 
1' ,/ ,10X,52 ('-')) 

WRITE (6,544) NCASE,P,R,NOBS,NVAR 

544 FORMAT (/,10X, 'CASE NO.=' ,I5.5X, 'P=' ,F7.4,5X, ’R=',F7.2,5X, 'N=' ,Ir,5 
IX, 'NO. OF VALUES=’ ,13,/) 

SUM=0 . 0D0 
SUM1=0 . 0D0 
DO 445 I=1,NVAR 
SUIYNSUM4FRQ (I) *DAT (I) 

445 SUM1=SUM1+FRQ (I) *DAT (I) **2 . 0D0 

XBAR=SUM/DFLQAT (NOBS ) 

SSQR= (SUM1-DFIOAT (NOBS) *XBAR**2 . 0D0) /DFIOAT (NOBS-1) 

IF (XBAR.GE.SSQR) GO TO 71 
P=XBAR/SSQR 

R=XBAR**2 . 0D0/ (SSQR-XBAR) 

GO TO 72 

71 R=5 . 0D0 
P=0.9D0 
GO TO 114 

72 Q=1 . ODO-P 
K=1 

7 TEST=Q** (K+l) / (DFLQAT (K+2) ) 

IF (TEST.LE.EPS) GO TO 8 
K=K+1 

GO TO 7 

8 M=K 
DETLAM+0 . 0D0 
DO 10 N=1,M 
RATICN=1.0D0 
DO 12 JJ=1 , N 

12 RATICW=RATION* DFIOAT ( JJ) / (DFIOAT ( JJ) +R) 

DEITAM=DETLAI4+0* *N*RATIOSI/DFLQAT (N+l ) 

10 CONTINUE 

L=1 

95 CALL SVS (FCN,TNOT) 

EFF1=-FCN (TNOT) 

RSTAR=R 
PSTAR=P 
SUM2=0 . 0D0 
DO 145 J=1,NVAR 

145 SUM2=SUM2+FRQ(J) (TNOT**DAT(J) 

TXBAR=SUM2/DFLQAT (NOBS ) 

CALL CQMPR(TNOT,TXBAR,XBAR,RSTAR) 



63 



P=R/ (R+XBAR) 

Q=1 . ODO-P 
K=1 

77 TEST=Q** (K+l) / (P*DFLQAT (K+2) ) 

IF (TEST.LE.EPS) GO TO 88 
K=K+1 
GO TO 77 
88 M=K 

DETLAM=0.0D0 
DO 51 N=1,M 
RATTON=l . ODO 
DO 52 JJ=1,N 

52 RATIOtt=RATION*DFLC!AT ( JJ) / (DFLOAT ( JJ) +R) 

51 DETLA^DETLAMfQ* *N*RATION/DFDOAT (JSH-1 ) 

EFF2=-FCN (TOOT) 

IF (L.NE.l) GO TO 141 
UNOT=-DLOG (DABS (TOOT) ) 

WRITE (6, 98) PSTAR, RSTAR, P , R, EFF1 , EFF2 , UNOT 
98 FORMAT (/, 11X, 57 ('&') ,/,15X, 'VALUES OBTAINED AT THE END OF FIRST IT 
1ERATION' ,/,12X, 'PTHIL' ,3X, 'RTHIL' ,3X, 'PSTAR' ,3X, 'RSTAR' ,3X, 'ETHIL' 
2,3X, 'ESTAR' ,3X, 'UNOT' ,/,12X,F5.3,2X,F7.3,2X,F5.3,2X,F7.3,2X,F5.3,3 
3X, F5.3,2X,F7. 4, /,11X,57 ('&'),/) 

141 IF (DABS (P-PSTAR) .GT. 0.001) GO TO 29 
IF (DABS (R-RSTAR) .GT. 0.005) GO TO 29 
IF (DABS (EFF1-EFF2) .GT. 0.01) GO TO 29 
UNCT=-DL0G (DABS (TOOT) ) 

WRITE (6,94) XBAR,SSQR,P,R,EFF2,UNCT,L 
94 FORMAT (15X, 'USING THE ABOVE DATA FOLLOWING INFORMATION IS OBTAINED 

1/ ,/ ,15X, 54 ( '-' ) ,/,10X, 'SPLE.MEAN' ,2X, 'SPLE.VAR' ,3X, 'EST.QF P',3X, ' 
2EST.OF R' ,6X, 'EFF' ,6X, 'OP. UNOT' ,3X, 'NO.CF ITR' ,/ ,10X,7 (9 ( '-' ) , 2X) , 
3/,10X,F9.3,2X,F9.3,2X,F9.3,2X,F9.3,2X,F9.3,2X,F9.4,5X,I3,//) 

GO TO 114 
29 L=LHL 

GO TO 95 
114 CONTINUE 

113 CONTINUE 

112 CONTINUE 

STOP 
END 
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SUBROUTINE GNRNB(N,0BS,P,R,ISEED2,ISEEDl) 
REAL* 8 P,R,OBS 

DIMENSION OBS(200) ,GAM(200) ,K(1) 
P1=SNGL(P) 

R1=SNGL(R) 

CALL GAMA (R1,ISEED1, GAM, N) 

DO 1 1=1, N 

RLAM=GAM(I) * (1.-P1)/P1 

CALL GGPOSH (RLAM, ISEED2 , 1 , K, IER) 

1 OBS (I) =DBLE (FLOAT (K (1) ) ) 

RETURN 

END 



SUBROUTINE CQMPR(T,TX,XB,RST) 

IMPLICIT REAL* 8 (A-H,0-Z) 

COMMON P,R,DETLAM 
R1=RST 

30 B=R1+XB* (1. 0D0-T) 

TY=DABS (Rl/B) **R1 
F=TY-TX 
FP1=TY 

FP2= (XB* (1 . 0D0-T) /B) +DLOG (DABS (Rl/B) ) 

FP=FP1*FP2 

R2=R1-F/FP 

IF (DABS (R2-R1) .LE.0.0001D0) GO TO 66 
R1=R2 
GO TO 30 
66 R=R2 

RETURN 

END 
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